Ground-state phase diagram and critical temperature of two component Bose gases 

with Rashba spin-orbit couphng 
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Ground-state phase diagram of two-component Bose gases with Rashba spin-orbit couphng is 
determined via a variational approach. A phase in which the fully polarized condensate occupies 
zero momentum is identified. This zero-momentum phase competes with the spin density wave 
phase when interspecies interaction is stronger than intraspecies interaction, and the former one is 
always the ground state for weak spin-orbit coupling. When the energies of these two phases are 
close, there is a phase separation between them. At finite temperature, such a zero-momentum 
condensation can be induced by a ferromagnetic phase transition in normal state. The spontaneous 
spin polarization removes the degeneracy of quasiparticles' energy minima, and consequently the 
modified density of state accommodates a Bose condensation to appear below a critical temperature. 

PACS numbers: 67.85.Fg, 03.75. Mn, 05.30.Jp 
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The experimental progress of spin-orbit (SO)-coupled 
quantum gases has generated great interest recently ^- 
|4]. Unlike the SO coupling extensively studied in solid 
matters, the synthetic SO coupling in dilute cold atomic 
gases can be realized in feasible configurations with a 
tunable coupling strength [5j 6^. Moreover, it provides 
a unique platform to study the rich physics of SO ef- 
fects in a bosonic system [7^12] . which have not been 
explored before. One particular interesting configuration 
is Rashba SO coupling [13l[T4j. In free space, since the 
minima of the single-particle spectrum have infinite de- 
generacy, the ground state of a Rashba SO-coupled Bose 
gas is in fact determined by the details of interactions. 
Two different Bose-Einstein condensation (BEC) phases, 
i.e., the plane-wave (PW) phase and the spin density 
wave (SDW) phase (also referred to as the stripe phase), 
have been identified in a previous study [5]. In a har- 
monic trap, exotic half-vortex states appear for suitable 
interaction parameters pTOJ. At finite temperature, al- 
though a condensate can not survive in a noninteracting 
system because of the constant density of state (DoS) 
in the low-energy limit, the BEC state is indeed ener- 
getically favored in the presence of repulsive interactions 
[T^ . In this sense, interactions do help the emergence of 
condensation. 

Previously, many theoretical works focus on the unique 
degeneracy of the single-particle spectrum. However, in 
a many-body system, the condensate wavefunction is not 
necessarily constructed by the single-particle states with 
minimal energy. In this work, we present a different 
ground-state phase diagram as summarized in Fig. 1. In 
addition to the PW phase and the SDW phase, a fully 
polarized phase in which the condensate occupies zero 
momentum (ZM) is identified, which is in strong con- 
trast to single-particle physics. At finite temperature, 
the emergence of ZM condensation can be induced by a 
ferromagnetic phase transition in the normal state. The 
spontaneous spin polarization generates an effective Zee- 
man field, which significantly modifies the spectrum and 
the DoS of quasiparticles. Consequently, a finite BEC 
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FIG. 1; (Color online) Ground-state phase diagram of Rashba 
SO-coupled Bose gases. The three BEC phases, i.e., PW 
phase, SDW phase, and ZM phase, have distinct magnetic 
properties as indicated by the spin polarization Af . There is 
a phase separation (PS) in the shadow region. 



temperature is retained. 

We start from a single-particle 
Rashba SO couphng {h = 1), 



Hamiltonian with 



HsO = 7, 1 {Px(Jx+PyCry), (1) 

2m m 



where ko is SO coupling strength, and a^^y are Pauli ma- 
trices. The eigenstates of this Hamiltonian construct two 
helicity branches with dispersion Ep-t = e^ ± k^p^/m. 
Here, Cp — p^ /{2m), and quantum number helicity a — zt 
denotes spin parallel or antiparallel to the in-plane mo- 
mentum pj_ . The minima of the energy spectrum are lo- 
cated at the ring in momentum space with p± = kg. For a 
noninteracting Bose gas, the ground state has infinite de- 
generacy, and the condensation wavefunction could be an 
arbitrary superposition of the single-particle states with 
the momenta on the degenerate ring. 

For a dilute system, interactions between atoms can be 



described by contact potentials, 

dr Qaa'i^l'^^lii^a'^a (2) 






with tAo- being the field operator and g^^ — g^^. For 
simplicity, we consider g^^ = g^^ = g m. this work [15] 
and assume all the interactions are repulsive. At the 
mean-field level, the interaction energy for a BEC state 
is given by 

fint = \j dv [{g + gnW{v)~5gnH'P{v)\ , (3) 

where 5g = g-fi - g, n ^ n-f + n^, and M = (ji^ - n^)/n 
are local density and local spin polarization respectively, 
and n = N /V is average density. 

Even without a serious calculation, a qualitative predi- 
cation can be made from Eq. (pi . Considering the 5g > Q 
case, interactions prefer either a fully polarized state or 
a SDW state to avoid the interspecies repulsion, and the 
former one has a lower interaction energy. On the other 
hand, kinetic energy is controlled by the parameter /cq. 
Therefore, when SO coupling is weak enough, kinetic en- 
ergy can be safely ignored, and the ground state is a 
fully polarized phase. This argument actually explains 
an important feature of the phase diagram in Fig. 1. 

T ground-state phase diagram can be determined 
quantitatively via a variational approach. For a state 
with uniform density and polarization, the condensate 
wavefunction is generally written as 



V 



cos 7 
' sin 7 e*' 



(4) 



where k is the condensation momentum, (j)\^ = arg(fc2: + 
iky) for fc 7^ and an arbitrary number for fc = 0, 
and the parameter 7 e [0, 7r/2] determines the spin po- 
larization. The variation energy is readily obtained as 
£ = N{tk + \9n- ^fcj_sin27-f jj^nsin^ 27). By mini- 
mizing this energy with respect to parameter 7 and mo- 
mentum k, we find that when Sgh < Aek^ , an unpolarized 
condensate with (k±,kz) = (fco,0) is preferred, and the 
wavefunction in Eq. Q corresponds to the single-particle 
energy minima. This state is called the PW phase, whose 
energy is given by, 



fpw ^N 



(m + 9)n - e/co 



(5) 



When 6gn > Aeko , the fully polarized condensate which 
occupies zero momentum is energetically favored, i.e., 
k = 0, M = cos 27 = ±1. The energy of this ZM phase 
is given by 



•fzM = ^Ngn. 



(6) 



We note that condensation at ZM is a pure interaction 
consequence, which is in strong contrast to the single- 
particle picture. In an equivalent interpretation, the 



spontaneous spin polarization generates an internal ef- 
fective Zeeman field /loff = SgnM/2, which modifies the 
energy minimum. As discussed later, a similar effect also 
appears in the normal state above BEC temperature. 

For a BEC state that breaks translational symmetry, 
condensate wavefunction can be written as (k ^ 0) 



V> 



C0S7 
-sin7e 



J0k 



COS7 
-sin7'e*'^-'' 



with |cp -t- |c'p = 1 and 7,7' £ [0,7r/2]. The ansatz of 
the above variational wavefunction is a minimal model to 
accommodate with spin and/or density spacial modula- 
tion. At the mean-field level, only kinetic energy depends 
on momentum k, and a straightforward variation yields 
fc^ = and k_i = ko{\c\'^ am2-f + |c'psin27'). Thus, to- 
tal energy £ can be expressed as a function of |cp, |c'p, 7, 
and 7'. Although an analytic solution to this variational 
problem is difficult, numerical energy minimization is 
straightforward by exhausting all configurations in pa- 
rameter space. We find that when a non uniform solution 
has a lower energy than the PW phase and the ZM phase, 
condensate wavefunction always takes a SDW form: 



<P 



^iri 



1 



^ik-r 



1 



-ik-r 



(7) 



where {k±,kz) = {ko,0), and the relative phase rj deter- 
mines the node position of the spin density modulation 
M(r) = cos(2k-r-)-?7). The SDW condensation described 
by Eq. ^ was first discussed in Ref. [S]. 
The energy of the SDW phase is given by 



^SDW = N 



35 + gn - 



-n - ekr, 



(8) 



By comparing this energy with Eqs. (pi) and ([6]), we ob- 
tain the ground-state phase diagram as shown in Fig. 1. 
When Sg < 0, the PW phase is energetically favored; 
when Sg > 0, the ground state is either in SDW phase 
or ZM phase. The competition between these two phases 
is controlled by SO coupling strength k^, and a first- 
order transition takes place at 6gn = 8efcp (indicated by 
the dashed-line in Fig. 1). When SO coupling is weak 
enough, the ZM phase is always the ground state (for 
Sg > 0) [Tg. 

Up to now, we assume the ground state is homoge- 
neous. Since the condition for the transition between 
the SDW phase and ZM phase is density dependent, 
a phase separation (PS) may appear. The boundaries 
of the PS region are determined by the balance condi- 
tions for chemical potential /i = d£/dN and pressure 
P = {iiN — £)/V of these two phases, 

MSDw("-i) = A^zm("2), PsBwini) = Pzu{n2)- (9) 
From the equation of state in Eqs. ^ and ((sl), we find 

(10) 



Til 



n2 



35 + .9n - 2^.9(3.9 + .9n)' 

2efcn 
^5(35 + m) - 2g 



(11) 
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FIG. 2: (Color online) Elementary excitations of the ZM 
phase for giifi = 2(7n = 20efc(,. The spectrum has a rota- 
tion symmetry in the Px-Py plane and a reflection symmetry 

Ep{6) ^Ep(-K-9). 



For the density h € (711,712), a mixture of the ZM phase 
and the SDW phase has a lower energy than in the ho- 
mogeneous case, and thus a phase separation takes place. 
The PS region is shown as the shadowed area in Fig. 1. 
When gn ^ e^^ , the PS region becomes very narrow, 
and the boundaries approach the first-order transition 
line Sgn = Sckg ■ 

Beyond the mean-field level, one can use Bogoliubov 
theory to take account of fluctuations. Previously, ele- 
mentary excitations of PW phase and SDW phase have 
been studied [H [T71 [TS]. Here, we focus on the ZM 
phase. 

Assuming a pure spin-up condensate, the Bogoliubov 
Hamiltonian for the ZM phase is given by 



/CzM - J2 [H^^P - i^^P 



-9n' 



(12) 



with *J, = ii^l^^, 



t 



/ ^P + 



gn 



,t 
fcp 



,'0-p,t'V'-p,i) and 



p±e 
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gn 




cp 
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gn 
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'^p^e-"*'- 



rP-Le' 



Sgn / 



and it reaches a minimum value when the phonon mode 
propagates in the Px-Py plane. In the /cq ~^ limit, 
Eq. (jTSJ) recovers the usual sound velocity of a single- 



component BEC. 

Now, we consider the normal state without a Bose con- 
densation. 

In a noninteracting system, the equation for total den- 
sity at temperature T is given by (fee — 1) 



p,a— ± 



p,aj 



)= deD{e)f{e), 



(14) 



where /(e) — [e^'^^'*-'/"^ — l] is Bose distribution func- 
tion, and D{e) = y ^ ^ 5{e — ep^a) is the DoS of single- 
particle. Due to the degeneracy of energy minima, D{e) 
is a constant in the low-energy limit [H] , and such a DoS 
prevents a condensate surviving at finite temperature. 



In the presence of interactions, Eq. (14) is still valid 



at the mean-field level with ep^± and D(e) replaced by 
the dispersion and the DoS of quasiparticles respectively. 
If Hartree-Fock energy only makes a constant shift in 
the quasiparticle spectrum, the DoS remains the same as 
in the noninteracting case, and the normal state is not 
kinetically forbidden at any finite temperature 12J. 

The situation is different once a spontaneous magneti- 
zation appears. The mean-field Hamiltonian for a normal 
state with a possible spin polarization M is given by 

P cr=tl 



+ (^P^^"^''^lJp,t + h.c. 



fP 

"^int! 



with £j'Jj^ = V{gn'i + gnj^ + g-fin-^nj^) and a the opposite 
spin to a. The quasiparticle energy spectrum ep.± = 
ep + {2g + g')n± yJ{g'nMY + ipj_ko/niy has a different 
structure from the noninteracting case for Af 7^ 0, where 
the interaction parameter g' = g^i/2 — g is defined. 
Spin polarization M should minimize the free energy 



J'n = /iA^ + r Y^ In (1 



p(M-£p.»)/T^ 






(15) 



p,Q — ± 



and the stationary point condition dJ-^/dM — can be 
explicitly written as 



By diagonalizing Eq. ( |12[ ), two branches of elementary 
excitations i?p,± are obtained. Some examples of -Ep,± 
along different directions are shown in Fig. [2J Within 
the ZM region of the phase diagram, we find the spectra 
are always positive. This result disagrees with those of 
previous works I19j , which claim the zero-momentum 
condensate always suffers a dynamical instability I20j. 

The higher excitation branch i?p,+ possesses an energy 
gap 5gn, which is just the energy cost to flip a spin to 
overcome the effective Zeeman field. The lower branch 

Ep _ is gapless with a phonon dispersion in the long wave which is equivalent to the self-consistency equation M = 
length limit, hmp^g^p,- = cp. The anisotropic sound T^Hp {Wp,-\i^T>.\) ~ i'i'p,ii'p,i))- For a given tempera- 



M 



-E 

V ^ 



ff-Af[/(£p,_)-/(£p, + )] 

^{g'nMY + (pifco/m)2 



(16) 



velocity c depends on direction angle 9 = arccos(pz/p), 



c{e) 



4£fco sin^ e 

{an - 9)n 



(13) 



ture, M and /i can be solved from Eqs. (14 1 and (16). 

For g'n < 2ekg, Eq. (16) only has a trivial solution 
M — ai any temperature, which means the normal 
state is always unpolarized. For g'n > 2efcQ, a finite spin 
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FIG. 3: (Color online) (a) Ferromagnetic transition temper- 
ature Tni and condensation temperature Tc as functions of 
'h}'''^a' . (b) Temperature dependence of chemical potential 
H for n^/'-^a' = 0.05 (D), 0.06(a), and 0.07(o) with fixed 
fco = 0.25n, and the dashed line is for the noninteracting 
case. The inset is the corresponding spin polarization M 
as a function of T. (cl)-(c4) Energy spectrum of quasipar- 
ticles Ep,- for the data labeled by • in panel (b), and the 
dashed lines indicate the value of /i — £min for each case. Here, 

12/3 



T^^2^[n/al)]'"/m. 



magnetization is energetically favored below a ferromag- 
netic transition temperature Tm, which is determined by 



1 



^s:,-5i[^'=-'-^<^-fl- 



V 



(17) 



with M = 0. In Fig. Is] (a), Tm is plotted as a func- 
tion of n^^'^a' for different /co/n, with a' = g'-m/iA-n:). 
One can see that Tm increases rapidly when g'n beyond 
the threshold 2tka- For small fcg/n, the phase transition 
takes place in the weak interacting regime {n^/^a' <C 1), 
where the validity of mean-field treatment can be jus- 
tified. We note that similar ferromagnetic phase tran- 
sitions also happen in the normal state of usual spin-1 
bosons [22] and spin-^ bosons [23] without SO coupling. 
Below Tfn, M increases very quickly as shown in the 
inset of Fig. [s] (b). The spin polarization generates 
an effective internal Zeeman field /icff = g'nM, which 
modifies both the spectrum and the DoS of quasiparti- 
cles. For M < 2ek^J{g'n), quasiparticle energy minima 
still construct a degenerate momentum ring (with a re- 



duced radius), and DoS remains a constant in low energy 
limit. However, when M > 2€ko/{g'n), i-e., hes > 2efco, 
zero momentum becomes the solo minimum of the quasi- 
particle spectrum, and the DoS vanishes in low energy 
limit 24\. Such a DoS accommodates a condensate to 
appear, when chemical potential reaches the bottom of 
quasi-particle spectrum e„iin [see Figs. [3tcl)-|3tc4)]. Be- 
low the condensation temperature Tc, the normal state 
does not exist. As shown in Fig. [sla), an experiment- 
achievable Tc can be realized for a reasonable interaction 
parameter n^'^a' . 

In contrast to the noninteracting case, the finite BEG 
temperature obtained here is due to the interaction ef- 
fect. We emphasize that the emergence of ZM conden- 
sation is induced by a ferromagnetic transition in the 
normal state, and the critical temperature T^ is natu- 
rally derived according to Einstein's original idea, which 
is quite different from the energetic scenario proposed in 
Ref [12 . On the other hand, for g'n < 26^^, the con- 
densation transition must be determined energetically 
(at the mean-field level), because normal state the is 
not kinetically forbidden at any finite temperature. In 
the T —>■ limit, we find the free energy of the un- 
polarized normal state for g'n < 26^^ is given by [25j . 
-Fn(T ^ 0) = N[i2g + 5n)"/4 - efeo]> which is in- 
deed higher than the energy of the BEG state given by 
Eqs. ([5|, (|6]), and ([8|. Therefore, a phase transition be- 
tween the normal state and the BEG state is expected at 
the temperature when their energy levels cross |12j . 

At intermediate temperature < T < Tc, a com- 
plete phase diagram is still unclear at present. In prin- 
ciple, the temperature dependence of spin polarization 
and condensate fraction can be studied via a Hartree- 
Fock-Bogoliubov approach. However, this problem be- 
comes complicated when the competition between the 
ZM phase and the SDW phase is considered. We leave 
this issue to future study. 

In summary, a spin-polarized ZM state is identified in 
the phase diagram of Rashba SO-coupled Bose gases, and 
the emergence of ZM condensation at finite temperature 
can be induced by a ferromagnetic transition in the nor- 
mal state. The phase transition and the critical temper- 
ature studied here may be useful to future experiments. 

Note Added. Recently, K. Riedl and coworkers studied 
the magnetic transition of the PW phase with g^^ < 
|26 , and a finite condensation temperature was also ob- 
tained. 
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